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Abstract 

^ ' We investigate a number of observables that are and will be instrumental in the exploration of tt pro- 
'— duction and decay at the Tevatron and the LHC. For this analysis we made a computer program that 
incorporates besides the NLO QCD corrections to tt production and decay also mixed weak-QCD cor- 
J> rections to the production amplitudes, and that allows for studies of correlated versus uncorrelated tt 
■ events. In this paper we analyze and compute observables mostly for dileptonic tt final states to next-to- 
, leading order in the strong and weak gauge couplings (NLOW), with selection cuts. We calculate charge 
CO I asymmetries of the top quark and of i = e, and we compare, where possible, with experimental results. 
CO ' We show that top-quark spin correlations affect a leptonic pair asymmetry (which has not been measured 
: so far) by ~ 7%. We determine several dileptonic angular correlations, which reflect tt spin correlations, 
_ namely the beam, off-diagonal and helicity correlation, and the opening angle distribution (defined in the 
^ ' t, t rest frames) when selection cuts are applied. Our NLOW predictions for the beam, off-diagonal, and 
•'H helicity correlation for the Tevatron agree with recent measurements by the CDF and DO experiments. 
rN ■ In addition we make predictions for estimators of these correlations as functions of M^f . These estimators 
' may prove useful for the Tevatron and also in the early rounds of LHC data- analyses, where the event 
numbers will not be abundant. Furthermore, we compute to NLOW in the gauge couplings two dilepton 
angular correlations that are defined in the laboratory frame, for correlated and uncorrelated tt events 
at the LHC. Finally, based on our SM results, we identify several observables that allow to search for 
non-standard parity- and CP-violating interactions, especially with future LHC data. 
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I. Introduction 



Exploration of the production and decay of top-quark pairs (tt) is and will be among the central 
physics topics at the Tevatron and at the Large Hadron Collider (LHC). On the theoretical 
side, the investigation of these processes has reached quite a sophisticated level. The next-to- 
leading order (NLO) QCD corrections to the ti cross sections and several top-quark distributions 
have been known for quite some time from the pioneering work of and [SI [6], which was 
extended by threshold resummation calculations [714T0]. The pair cross sections were recently 
updated by P^HTT] . The NLO QCD corrections including the full dependence of the matrix 
elements on the top-quark spin degrees of freedom were determined in [TTHIl]. The mixed weak- 
interaction QCD corrections [T8]-[2l]and the photonic corrections |25] to tt production and the 
order nonfactorizable QCD corrections |26] are also known. The computations of the NLO 
QCD corrections to ti + jet, ti + bb, and ti + 2 jets, which are important background processes in 
the search for (the) Higgs boson(s), were reported in [271EH], in [29143T] . and in [32] respectively. 
Differential distributions for polarized semileptonic and non-leptonic decays of the top quark were 
determined order by [33] and [31] , respectively, while its total width is known to order [35|l36] . 
A number of observables that are instrumental in exploring the ti production dynamics were inves- 
tigated in detaifl These include the VT-boson helicity fractions in t — )■ bW at NLO in the gauge 
couplings [10] and to QCD [H], and the top-quark charge asymmetry at the Tevatron in the 
context of standard model (SM) interactions [21I11I2ZII12H1S] • This asymmetry was measured by the 
DO and CDF experiments [TH - HO] . Top-quark spin effects, in particular ti spin correlations were 
analyzed in [50] - [6l] at LO and in [T3|[Ti] at NLO QCD. Measurements of ti spin correlations at the 
Tevatron were reported in [651467] . Quite recently there has been an increased activity on theoret- 
ical and phenomenological issues, including the first determination of the running mass frit of the 
top quark [68j from c^^^, analytic phase-space integration of the NLO QCD parton cross sections 
99 1 1'll 91 ~^ [HH], threshold expansion improvements of the partonic cross sections [701 - I73] 
and of the pair-invariant mass distribution [71], an analysis of quasi-bound-state effects in au at 
the production threshold [75|[76]. determination of various building blocks for the computation of 
att at NNLO QCD [771483] . phenomenological analyses of boosted top-quark events [8ll|85] and of 
the use of the LHC ti cross section for calibrating parton distribution functions [HS]- The NLO 
QCD Monte-Carlo generators [SZIISB], [SS], and [90j are important tools for the simulation of ti 
production and decay including parton showering [871|90]. Recently, results were presented based 
on a code which incorporates NLO QCD corrections to ti production and semileptonic top-quark 
decay with ti spin correlations [91] . 

In this paper we analyze a number of distributions and correlations, mostly for dileptonic ti par- 
tonic final states, at next-to-leading order in the strong and weak couplings. The spin-dependent 
terms are included to this order of the perturbation expansion in the production and decay matrix 
elements for the ti intermediate state. For our analysis we made a computer code with input 
based on previous work [T? ] [T1 |[TU]4^ IM] where the authors of this paper were involved. As our 
approach is based on spin density matrices being independent of specific spin reference axes, we 
can switch off the t- and/or t-spin dependent terms in the matrix elements, in particular the ti 

^For more complete overviews, of., e.g., [57H5^ . 
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spin correlations, also at NLO in the gauge couplings. This seems to us quite useful for perform- 
ing comparative studies of final state distributions for correlated and uncorrelated tt events. This 
paper addresses the following issues: i) A study of the size of the weak interaction corrections 
with respect to NLO QCD corrections for several distributions and spin correlations. So far, only 
comparisons with respect to LO QCD results were made in the above-mentioned literature, ii) An 
investigation of the top-quark charge asymmetry and of two leptonic charge asymmetries for 
correlated and uncorrelated dileptonic ti events at the Tevatron if acceptance cuts are applied. To 
our knowledge, this has not been analyzed so far to NLO in the gauge couplings, iii) A detailed 
analysis of dileptonic angular correlations induced by ti spin correlations. 

In Sect. [Ill we describe the ingredients of our analysis and list our input parameters. Sect. UTTl 
contains results for several distributions and spin correlations at the level of the ti intermediate 
states. The top-quark charge asymmetry at the Tevatron will be considered in some detail and 
a comparison with the results of |15] will be made. In Sect. IIVI we give our results for several 
distributions and angular correlations for dileptonic final states with phase-space cuts at NLO in 
the strong and weak couplings. Two distributions will be computed for semileptonic final states. 
We analyze the top-quark charge asymmetry and a related pair asymmetry at the Tevatron and 
two related leptonic charge asymmetries for correlated and uncorrelated dileptonic ti events, and 
we compare the top-quark asymmetries with experimental results |Tn - H9] . Our predictions of 
the leptonic asymmetries have to await confrontation with measurements. Then we investigate 
a number of dileptonic angular correlations, for the Tevatron and the LHC at ^/s = 10 and 14 
TeV, at NLO in the gauge couplings that are (in)sensitive to ti spin correlations. Moreover, we 
compare with the recent Tevatron measurements [651467] . Based on our SM results we identify, 
in addition, several observables that allow for searches for non-standard parity- and CP-violating 
interactions, especially with future LHC data. Sect. |V] contains a summary and outlook. 



II. Theoretical set-up 

We analyze hadronic top quark pair production and decay, at NLO in the QCD and weak couplings, 
in the on-shell approximation for the t, i quarks, taking their spin degrees of freedom in the 
production and decay stage into account. To be specific we consider, at the parton level, the 
following reactions: 

gg,qq ^ bb + Af, (ILl) 
gg,qq^bb + Af + g, (11.2) 

g + q{q) ^bb + 4f + q{q), (ILS) 

where f = q,i,Vf. As the top quark is a narrow resonance, we employ the narrow width ap- 
proximation Tt/rrit — )■ and incorporate the factorizable QCD and mixed weak-QCD corrections, 
which form gauge-invariant sets. In this approximation the squared matrix element |A^p of the 
respective reaction is of the form 

2 1 r -n 1 TT^ 

' ' " ^*2l^^'^^ ^^^^^ " iVi — ---2^p„,„i?„„,,;3^,p^,^. (II.4) 
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Here R, p, p are the production and t and t decay density matrices for ( jll.ip - (]II.3|) . The subscripts 



in Eq. ( 1II.4P denote the top and antitop spin indices, and the factor iVj arises from averaging over 



the colors and spins of the initial partons. 

As already mentioned above, in our spin density matrix approach the t, i spin-dependent terms, in 
particular tt spin correlations can be switched off also at NLO in the gauge couplings. In addition, 
top-spin effects can be analyzed for arbitrary spin bases at NLO, in particular at the level of the 
tt intermediate states. This is quite useful for studies of non-standard interactions that change 
the SM-induced top-spin correlations. 

As to the density matrices R of the ti production reactions qq — )■ tt{g), gg — )■ ti{g), gq — )■ tiq, and 
gq — >■ tiq: Besides taking into account the 0{a^g) QCD corrections [11], we incorporate also the 
weak interaction corrections of order a^a, agOt-, o? ^ and UgCp' to the spin density matrices of these 
processes as computed in [T9H2T]. 

We use the decay density matrices of the main semileptonic and non-leptonic SM decay modes of 
polarized (anti)top quarks, t — )■ hW^ — )■ ht^V({g\ bqq'{g) (where qq' = ud,cs) to 0{as) in fully 
differential form. Specific distributions were given in [5^ . 

In the numerical implementation of (1II.4P we use the expanded form of the right-hand side to NLO 
in the gauge couplings. We expand also the top width Fj in the denominator of (III.4|) to 0{as). 
In Sections IIIII and IIVI below, the acronyms LO and NLO refer to the 0{a1) and 0{al) matrix 
elements in QCD, while NLOW refers to the additional inclusion of the above-mentioned weak 
corrections. 

The infrared (i.e. soft and coUinear) singularities in the production density matrices are taken 
care of by a phase-space slicing method, as worked out explicitly in ^14]. For semileptonic top- 
quark decay we use a slicing procedure, too, while for non-leptonic decay a hybrid scheme is 
employed (cf. [21]). We checked that the results of Sections IIIII and HVl are independent of the 
slicing parameter Xmm for x^m ^ a few xlO~^. The results of these sections were obtained with 
r = 10"^ 

With these building blocks we developed a numerical program for computing observables to fixed 
order in the gauge couplings as described above, at the level of the leptonic and partonic final states 
that correspond to the dileptonic, lepton plus jets, and all-jets final states. As we are primarily 
interested in modeling "top as a signal", we investigate in this paper only distributions mostly 
for dileptonic and a few for semileptonic final states which have the highest signal-to-background 
ratios. 

pp, pp ^ tt + X ^ £+i jb n + X, (II.5) 
pp, pp -> ti+X -^t jb hji 32 + X, r ]i jb ji 32 + X, (II.6) 

where i = e, p, r, and jb (ji,2) denote (non) b jets. 

Our input parameters are as follows. Throughout this paper the top-quark mass rrit is defined in 
the on-shell scheme. Adopting the common, albeit not rigorously justifiable practice of identifying 
the experimentally determined top mass with the on-shell mass, we use for rrit the recent Tevatron 
average rrit = 173.1 ± 1.3 GeV [92j. Our main results, the charge asymmetries and the top-spin 
induced angular distributions and correlations given in Sect. IIVI are, in fact, not very sensitive to 
variations Arrit of the top mass by a few GeV. The QCD coupling is defined, as usual, in the MS 
scheme. We employ CTEQ6L1 and CTEQ6.6M parton distribution functions (PDF) for LO and 
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NLO computations [86j, which correspond to as{mz) = 0.130 and as{mz) = 0.118, respectively. 
The evolution of as = as,5 to some higher scale hr is made within 5-flavor QCD. In the parton 
matrix elements, which were computed in p31 within 6-flavor QCD, the MS coupling as,6(/^i?) is 
replaced by (ys,5{f^R) according to asfi/n = as^^/ir — |(as_5/7r)^ ln(m^//i|,) + 0(0:^5) [93]. For the 
top width we use Ff = 1.3882 — as{mt) x 0.8076 GeV. We use an on-shell mass rrih = 4.8 GeV 
for the b quark, while the masses of the lighter quarks and of the leptons are put equal to zero. 
In the evaluation of the weak corrections we use the QED coupling a{mz)=0.008 and for the W- 
and Z-boson masses and the ly-boson width F^y the measured values given in [94]. We put the 
CKM matrix elements that appear in the CKM-unsuppressed semi- and non-leptonic top-quark 
decay amplitudes equal to one, |Vt;,| = \Vqq>\ = 1. The Higgs-boson mass is put equal to ttih = 120 
GeV. The dependence of the results of Sections UTTl and HVl on rriH < 200 GeV is insignificant. 

III. Results for ttX: 

In this section we give results, at LO, NLO, and NLOW in the gauge couplings, for several 
distributions and correlations at the level of on-shell ti intermediate states, 

pp,pp ^ ti + X , (III.l) 

both for the Tevatron {^/s = 1.96 TeV) and for the LHC at center-of-mass energies ^/s = 10 and 
14 TeV. No phase space cuts are applied in this section. The factorization and renormalization 
scales are put equal, fip = fi^ = ^ and we choose /i = mt/2, m^, and 2mt. 

As was shown in [20], the weak- interaction contributions to the total cross section au at the 
Tevatron and at the LHC are marginal as compared with the NLO QCD corrections. On the 
other hand, for a number of distributions, for instance, the transverse momentum distribution 
of the (anti)top quark and the ti invariant mass distribution, the weak-interaction corrections 
are of potential importance at large energies due to large Sudakov logarithms. For the pt [Mu) 
distribution at the LHC these corrections amount to about —10% (—6%) of the LO results for 
~ 1 TeV {Mu ~ 2 TeV), and the respective ratios daweak/daio grow in magnitude for larger 
Pt and M^f [T91423] . Here we compare the weak and NLO QCD corrections for tt production, to 
wit, for the px and distributions and for several ti spin correlation observables at the LHC, 
and for the charge asymmetry at the Tevatron. 

Fig. [1] shows, for the LHC at ^/s = 14 TeV, the (un)normalized ti invariant-mass distribution 
at LO, NLO, and NLOW for /i = rrit. The analogous plots are displayed in Fig. [2] for the 
(un) normalized pt distribution of the top quark. In Fig. [3] the ratios of these distributions eval- 
uated at NLOW and NLO are plotted for yu = rrit. The weak-interaction corrections to the 
distribution are negative except close to 2mt and become larger in magnitude than ~ 2% with 
respect to the NLO QCD corrections for M^f > 1.2 TeV. For the pt distribution the respective 
ratio is also smaller than one, except in the extreme forward- and backward region [21], and the 
weak-interaction corrections grow to several percent in magnitude beyond pt ~ 500 GeV (not 
shown in Fig. [3]). A study of the scale dependence of the normalized pt and distributions in 
the range Pt ^ 1 TeV, Mu ^1-2 TeV shows that for mt/2 < fi < 2mt the scale variations are of 
the same order of magnitude as the weak-interaction contributions to the NLO QCD distributions. 
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Figure 1: The unnormalized (left panel) and normalized (right panel) tt invariant- mass distribution 
for the LHC (^i = 14 TeV) at LO (solid), NLO (dashed), and NLOW (dotted) for // = rut. 




Figure 2: The unnormalized (left) and normalized (right) pr distribution of the top quark for the 
LHC {y/s = U TeV) at LO (solid), NLO (dashed), and NLOW (dotted) for = m*. 
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Figure 3: The ratio of the Mu distribution evaluated at NLOW and NLO with fi = nit (left 
panel) and the corresponding ratio of the px distribution (right panel) for the LHC (14 TeV). 

Next we consider the asymmetry with respect to interchange of the the t and t charges, which 
is known to be relatively large at the Tevatron. This asymmetry is generated at NLO QCD by 
the interference of even and odd terms under t o t in the amplitudes for qq annihilation and, 
likewise, for gq and gq fusion [2|IH Il2] - H6] . (The contribution of the latter two processes to the 
charge asymmetry is one order of magnitude below the qq contribution.) The weak interaction 
corrections to the qq and flavor excitation processes contribute, too. The conventionally defined 
differential and integrated charge asymmetry, A{y) and A, are: 

where y denotes the rapidities yt, yi of the t and t quark defined in the laboratory frame, and 
N{y) = datt/dy. 

As the pp state at the Tevatron is a CP eigenstate in the laboratory frame, CP invariance of 
the SM correction^ implies that Ni{yi) = Nt{—yt), which in turn implies that A is equal to the 
forward-backward asymmetry of the top quark: 

J^Ntiy)- J^Ntiy) 

y>0 y<0 



^The CP-violating terms in the matrix elements for pp tt final state induced by the non-zero SM KM 
phase Skm in higher orders of perturbation theory are of course numerically irrelevant. 
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In the case of the Tevatron, CP invariance tells us also that A 



FB 



-A 



FB- 



Non-standard CP- 
violating interactions, if existent, may invalidate this relation by a small amount. As to ti produced 
by gg fusion: the Bose symmetry of the gg state precludes a contribution from this sub-process to 
the integrated asymmetries A, A^p^ - irrespective of whether or not the production density matrix 
Rgg contains P- and/or CP- violating pieces. 



The DO and CDF experiments [47|,38] and [15] considered the following observable which also 
reflects the charge asymmetry in tt production at the Tevatron: 



J N{Ay > 0) - / N{Ay < 0) 
/ N{Ay > 0) + jN{Ay < 0) 



(III.4) 



where Ay = yt — yi- This asymmetry is, for kinematical reasons, larger than A. 
It should be recalled that the numerators of the above asymmetries start at order in the 
QCD perturbation expansion. Thus, the above inclusive asymmetries are of order as- They were 
computed to this order, including weak corrections, by [ISlllS]- The order QCD corrections 
are not known. Also within our framework, the determination of these asymmetries is a LO 
calculation. However, we compute here the numerators of A and A** by taking into account the 
NLO respectively NLOW ti parton matrix element^ both from gg, gq{q), and gg initial states. 
(In this way, the computation of these asymmetries provides also a check of our code.) We adopt 
the following procedure. We evaluate the numerators at NLO and NLOW, respectively, with 
NLO PDF and the denominators at LO with LO PDF, using however the same value of in 
both numerator and denominator. (Evaluating the denominator including the NLO corrections 
would be in conflict with the perturabtion expansion.) We label our results with the acronyms 
NLO' and NLOW', respectively. 

We have plotted in Fig. Hlthe differential charge asymmetry ( ]III.2|) for the Tevatron at NLO' and 
NLOW' for a fixed scale, /i = m^. Tabled] contains our results for the integrated charge asymmetry 
A and the pair asymmetry A** for the for three different scales. 



Table 1: Results for the charge asymmetry A and the pair asymmetry A** at NLO' and NLOW' 
for fi = mt/2,mt, and 2mj. The acronyms NLO' and NLOW' are explained in the text. 





Tevatron 




mt/2 


mt 


2mt 


A (NLO') 


0.054 


0.049 


0.045 


A (NLOW) 


0.056 


0.051 


0.048 


A" (NLO') 


0.084 


0.076 


0.071 


A" (NLOW) 


0.087 


0.080 


0.075 



■^In the calculation of [43l|45] just the relevant ©(a^) terms that generate the asymmetries were taken into 
account in the numerator. That is why they used LO PDF also in the numerator. 
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Figure 4: The differential charge asymmetry A{y), eq. (lIII.2p . for the Tevatron as a function of 
the (anti)top-quark rapidity y, for /i = nit. The solid and dotted lines correspond to NLO' and 
NLOW, respectively. 

Our results for A and A** given in Table [1] agree with the results A = 0.051(6) and A** = 0.078(9) 
of [45]. The following remarks are, however, in order. We use a different set of PDF and a slightly 
different value of m^. The contribution of the mixed electroweak-QCD corrections to the charge 
asymmetries were determined in [13] for the qq initiated production process which, as found by 
these authors, increases the QCD asymmetry A by a factor 1.09. We have taken into account 
in addition the weak-interaction contributions to gq (g) — ?■ tiq (q) and a specific contribution to 
bb — )■ ti [21j. In Born approximation, the amplitude of the latter reaction consists of t-channel W 

boson exchange bb-^ti and s-channel photon, Z boson, and gluon exchanges WP-^ti. The terms 
in \M.{bb — )• tt)\^Q that involve the electroweak interactions are of order and a^a, respectively. 
They contain also charge-asymmetric pieces which contribute to A and A** with opposite sign 
as compared to the contributions from the order a^a terms in the qq — )■ ti matrix element. At 
present the size of these bb induced contributions can not be determined precisely in view of the 
rather poorly known 6-quark PDF. In any case, these additional terms diminish the overall size 
of the weak-interaction contributions, as the numbers of Table [1] show. 

The spread of the asymmetries in Tabled] due to scale variations certainly underestimates the true 
uncertainties of these predictions of the inclusive asymmetric^. A more realistic estimate of the 
theory uncertainty, namely ~ 30%, is provided by the calculation of A in [16] where NLL QCD 
threshold resummation corrections were included. Comparisons with the experimental results 
from the Tevatron will be made in the next section. 

^For ti -f jet final states the charge asymmetry was calculated to NLO QCD in |27) . 
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At the LHC, the initial pp state is an eigenstate of parity. Thus, A^p^ = A^p^ = in the laboratory 
frame, as long as only parity-invariant interactions - more general, only parity-even terms in the 
scattering operator - are taken into account. The parity-violating terms of the weak corrections 
appear in fact only in the t- and/or t-spin dependent terms of the partonic production density 
matrices and do, therefore, not contribute when making predictions for top quarks summed over 
their spins. As a consequence, the differential charge asymmetry A{y) at the LHC induced by 
the SM interactions must be symmetric with respect to ?/ = 0. Thus, in the absence of cuts, the 
integrated charge asymmetry A = 0. We have checked that our code reproduces this conclusion. 
Nevertheless, there are non-zero effects in suitably defined distributions |43[l45]. 



Next we consider top-spin induced effects. The polarization of the t or t quarks in the hadronically 
produced ti sample due to weak interactions or QCD absorptive parts is small [20l[58]- On the 
other hand, the QCD dynamics leads to correlations of the t, i spins already at tree level which are 
characteristic of the production mode. One may investigate ti spin correlations by the following 
expectation values of spin operators O: 



Oab = {^{Sfk){S-fh)) and api,= (^Si-Sf). (III.5) 



Here Sj, Sf are the top and antitop spin operators, and the unit vectors a and h are arbitrary 
reference directions, which may be interpreted as spin axes. As usual, the normalization of the 
averages is such that (1) = 1. We use the following spin bases [H] in the zero-momentum frame 
(ZMF) of the ti pair: 

a = — b = k, (helicity basis), (HI. 6) 

a = = p, (beam basis), (HI-''') 
a = = d, (off-diagonal basis). (HI. 8) 

Here k denotes the direction of flight of the top quark in the tt-ZMF and p is the direction of one 
of the colliding hadrons in that frame. The unit vector d, which we use in the convention of [T4] . 
defines the so-called off-diagonal basis [60]. The correlations Oab (C'spin) are equivalent to (a sum 
of) double spin asymmetries. At the Tevatron, the strength of the spin correlations induced by 
the SM interactions is largest with respect to the beam and off-diagonal basis. For the LHC the 
helicity correlation and Ogpin are good choice^. 

We consider the variation of these correlations with Mu, which we define by 



dO 



\j daO 5W{pt+PtY - Mu) . (III.9) 



dMu a 

We normalize to the integrated cross section. 

The distributions are displayed in Fig. [S] to LO, NLO, and NLOW for the beam and off-diagonal 
basis at the Tevatron. For the helicity correlation and for Ogpin the analogous plots are given in 
Fig. [6] for the LHC (14 TeV). The latter distributions pass through zerqj near Mu ~ 800 GeV. 



^For gg — ^ tt an optimal axis was constructed in |61j . It was reconsidered recently in [64) . Cf. the next section. 
^Their detailed behaviour can be read off from plots of scaling functions given in jl4) . 
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These figures show the known feature that the ti spin correlations (1111.50 at the Tevatron and 
LHC receive most of their contributions from tt pairs with rather low Mtj. In Figs. [7] and [8] the 
respective distributions are given for fi = mt/2, rrit, and 2mt. 

In order to quantify the size of the weak corrections we have plotted in Fig. [Hlthe ratios of dO/dMti 
evaluated at NLOW and NLO for Cgpin and Chei at the LHC (14 TeV). This figure shows that 
the weak-interaction corrections to Cspm are negative and < 4% with respect to NLO QCD for 
Mti < 800 GeV, while these corrections are smaller for (9hei- In the latter case, the increase of 
the magnitude of the ratio in the vicinity of Mu ~ 800 GeV is just a consequence of the NLO 
correlation passing through zero and changing sign afterwards. 

In Table E] the values of these variables, integrated over the Mti spectrum, are given at NLOW for 
three scales. The spin correlations will be analyzed further for dileptonic final states in the next 
section. There we will also compare the numbers of Table [2] with the results of [H] . 

Table 2: Results for several spin correlation variables at NLOW at the level of ti final states (no 
cuts) for the Tevatron and LHC. At the LHC the SM-induced ti spin correlations are very small 
with respect to the beam and off-diagonal bases and are therefore not given. Their SM values 
were given in p^[T^. 





Tevatron 


LHC (10 TeV) 


LHC (14 TeV) 




mt/2 


rrit 


2mt 


mt/2 


TUt 


2mt 


mt/2 


m 


2m 


^bcam 


0.777 


0.791 


0.804 














0off 


0.783 


0.798 


0.810 
















0.216 


0.218 


0.220 


-0.236 


-0.233 


-0.229 


-0.239 


-0.236 


-0.234 


Chel 


-0.358 


-0.368 


-0.376 


0.326 


0.326 


0.325 


0.327 


0.328 


0.331 
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Figure 5: The distributions dO^^eaxa/ dMft (left) and dOos/dMu (right) for the Tevatron at LO 
(sohd), NLO (dashed), and NLOW (dotted) for n^rut. 
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Figure 6: The distributions dO^pin/dMu (left) and dOh^i/dMu (right) for the LHC {^/s = 14 TeV) 
at LO (sohd), NLO (dashed), and NLOW (dotted) for = rrit. 



11 



10 r 




J I 1 1 1 I 1 1 1 I 1 1 1 ja L I 1 1 1 I 1 1 1 I 1 1 1 1 

400 600 800 1000 400 600 800 1000 



M„ (GeV) M„ (GeV) 

Figure 7: The distributions dO^^^^/dMu (left) and dOos/dMu (right) for the Tevatron at NLOW 
for II — 171^/2 (dashed), /i — rrit (sohd), and /i — 2mt (dotted). 




Figure 8: The distributions dO^pin/dMu (left) and dOhd/dMu (right) for the LHC {^/s = 14 TeV) 
at NLOW for /j, = mt/2 (dashed), n = rrit (solid), and /j, = 2mt (dotted). 
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Figure 9: The ratios of dO/dMu evaluated at NLOW and NLO with fi = rrit for Ogpin (left panel) 
and Chei (right panel) for the LHC (14 TeV). 

IV. Results for dilepton and semileptonic final states 

In this section we give results for several distributions, correlations, and asymmetries for dileptonic 
final states ( 111.61) . and for a few observables for semileptonic final states ( 111.51) . at NLO in the strong 
and weak couplings (NLOW). The radiative corrections were implemented in our computer code 
as described in Section [Tll 

We use the following cuts for the dileptonic final states {i = e, /i, /Ey denotes the missing transverse 
energy, and r] is the pseudorapidity) : 

Tevatron : > 20 GeV, \rie\ < 2.0, > 20 GeV, \r]j\ < 2.0, > 25 GeV, 



LHC : > 20 GeV, \r]e\ < 2.5, > 20 GeV, \r]j\ < 2.4, > 40 GeV. (IV.l) 



For the semileptonic tt final states we use 

Tevatron : > 20 GeV, \rie\ < 2.0, p^^ > 20 GeV, \rij\ < 2.0, > 20 GeV, 

LHC : p^T>20 GeV, \r]e\ < 2.5, ^ > 20 GeV, \7]j\ < 2.4, > 20 GeV. (IV.2) 

We put again fiji = fip = fi and evaluate the distributions for fi = mt/2, rrit, and 2mt. We 
calculated the observables below inclusively. We checked for some dilepton final-state distributions 
that the results do not change when using instead the kj_ jet algorithm [95J. This is to say that 
we checked an inclusive calculation against d-^vd'T^i jhii (NLO) and (^^i^e^^i^e jbjij (LO), where j 
denotes a gluon or light quark jet. 

Fig. [To] shows, for the dilepton channels at the Tevatron and LHC (10 TeV), the p^ distribution 
of the charged lepton for three scales /i. We have normalized to the respective cross section 
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Figure 10: The normalized transverse momentum distribution of d.'^ in the id channels at NLOW, 
for the Tevatron (left) and the LHC at i/i = 10 TeV (right). The solid, dashed, and dotted lines 
correspond to /i = mt,mj/2, 2mi, respectively. 

a{ti — 7- X) in order to reduce the scale uncertainties. Nevertheless, as is well-known, the scale 
variations become significant at large pf^. Fig. [TT] shows the analogous plots for the normalized ti 
invariant mass distribution for the semileptonic final states. The only purpose of showing these 
plots is to demonstrate that our code produces also established fixed order NLO QCD results. 
The Mtt distribution is, especially for the semileptonic channels, an important observable in the 
search for resonances that (strongly) couple to ti - for investigations within various models, cf., 
for instance, [9611101] . 



IV. 1. Charge asymmetries 



Next we compute, for the dilepton final states, the charge asymmetry A ( ]III.2p . and the pair 
asymmetry A** f lIII.4p with the above cuts. The rapidity distributions of the charged leptons. 



da 

dye+ 



da 

dVi- 



(IV.3) 



are defined with polar angles of £^ and the rapidities y£± in the laboratory frame. They are shown 
in Fig. [12] (left) for the Tevatron, computed at NLOW with the cuts (lIV.ip . In this figure we have 
actually divided Ni± by the dileptonic cross section at NLOW. 
One may define a differential leptonic charge asymmetry by 



Ay) 



N,+ {y)-N,-{y) 
Ne+iy) + Ne-{y) 



(IV.4) 
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Figure 11: The normalized ti invariant mass distribution for semileptonic final states at NLOW, 
for the Tevatron (left) and the LHC at ^/s = 10 TeV (right). The solid, dashed, and dotted lines 
correspond to fi = mt,mt/2, 2mt, respectively. 



Table 3: Results for the charge asymmetry A and the pair asymmetry A** of top quarks, and 
the leptonic charge asymmetry and the pair asymmetry A^^, for dileptonic final states at the 
Tevatron, with the cuts flIV.2|) . at NLO' and NLOW for fi = mj/2, rrit, and 2'mt. The acronyms 
NLO' and NLOW are explained in Section IIIII The acronym ti uncorrelated refers to the case 
when all t- and t-spin dependent terms in the spin density matrices are switched off at NLOW. 





Tevatron {{tt correlated)) 


Tevatron {tt uncorrelated)) 




mt/2 


mt 


2mt 


mt/2 


mt 


2mt 


A (NLO') 


0.053 


0.048 


0.044 


0.053 


0.047 


0.043 


A (NLOW) 


0.054 


0.049 


0.046 


0.054 


0.049 


0.046 


A" (NLO') 


0.074 


0.068 


0.062 


0.075 


0.067 


0.061 


A" (NLOW) 


0.078 


0.071 


0.066 


0.077 


0.070 


0.065 


A' (NLO') 


0.038 


0.033 


0.031 


0.037 


0.033 


0.030 


A" (NLOW) 


0.039 


0.034 


0.032 


0.038 


0.035 


0.032 


A'" (NLO') 


0.047 


0.042 


0.038 


0.050 


0.045 


0.041 


A'' (NLOW) 


0.048 


0.044 


0.040 


0.052 


0.047 


0.043 
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where y denotes the rapidities yi+ and yi-. This distribution is of course not identical to the 
differential top-quark charge asymmetry A{y). It is shown in Fig. [12] (right) for the Tevatron, 
computed at NLOW (i.e., by the procedure described in Sect. IIII|) with the cuts (lIV.ip . 
The associated leptonic charge asymmetry is 

/ N,^{y) - J N,-{y) 

J N,,{y) + J N,-{y) ' ^ ^"'^ 

y>0 y>0 

If CP invariance holds, then Nf+{y^+) = N^-{—yi-) and is equal to the leptonic forward- 
backward asymmetry, A^ = A^p^ = —Ap^, where 

/ Ne4y) - I Ne^iy) 
Af^ = . (IV.6) 

y>0 y<0 

In analogy to the tt pair asymmetry A** one may consider a leptonic pair asymmetry 

^ jNiAye>0)-jNiAy,<0) 

jN{Ay,>0) + jNiAye<0)' ^ ' ' 

where Ayi = yi+ — yi-. The observables (lIV.6p . (]IV.7p and related, more exclusive asymmetries 
were analyzed in [44J for dileptonic and semileptonic final states at the Tevatron with an approx- 
imate NLO QCD procedure (for details,see [11]). Here we compute (IIV.5I) and (1IV.7I) at NLOW 
with t, i spin effects included and, for comparison, also for uncorrelated tt events. For the evalua- 
tion of A, y4**, A^{y), A^, and A^^ we use the procedure described in the previous section. Results 
are given in Table [3l at NLO' and NLOW, for the dilepton final states at the Tevatron. 
The distributions fllV.Sp on which A^ is based are displayed in Fig. [12] (left). The figure shows that 
the CP-symmetry relation N(+{y) = N^-{—y) is satisfied. The figure displays also the consequence 
of the qq initiated SM top-quark charge asymmetry which tells us that the t {t) is preferentially 
emitted into the (anti)proton direction. The numbers in Table [3] show that the charge asymme- 
tries A and A** become smaller by several percent as compared to the values given in Table [1] 
due to the acceptance cuts. Kinematic arguments suggest that A^ has the same sign than A but 
is smaller in magnitude, while the SM value of A^^ is larger than A^, in analogy to A** versus 
A. This is is corrobated by the the numbers in Table [3] We give in this table also the values 
of these asymmetries when all the spin-dependent terms, in particular the ti spin correlations 
in the matrix elements flll.4p are switched off at next-to-leading order in the gauge couplings. 
This corresponds to the production of uncorrelated ti pairs, followed by their decays which are 
spherical in the t and t rest frames. One sees that the values of the inclusive asymmetries A and 
A^ remain essentially unchanged, which is not surprising. As to A**, one may expect an effect 
of spin correlations on this ordered pair asymmetry in the presence of acceptance cuts. Yet the 
numbers for the correlated case are larger than those for uncorrelated events by only about 1%. 
An inspection of the top-charge asymmetric terms in the partonic matrix elements reveals that the 
spin-dependent terms make indeed only a small contribution to this rather inclusive observable. 
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However, as the numbers for the leptonic pair asymmetry A show, there is a moderate effect of 
tt spin correlations on this more exclusive observable. The difference between our predictions for 
correlated and uncorrelated ti events amounts to about 7%. 

In view of theoretical uncertainty of 6A ~ 30% that is implied by the analysis of jlH] at the 
level of tt intermediate states, one may assign an uncertainty of this order of magnitude also to 
the numbers of Table [31 The investigation of more exclusive asymmetries at NLOW, e.g. those 
proposed in [H], including semileptonic tt events and the corresponding asymmetries at the LHC 
is beyond the scope of this paper and has to be left for a future investigation. 
The DO experiment at the Tevatron obtained, for lepton plus jets final states, a pair asymmetry 
A** = 0.12 ± 0.08 ± 0.01 |17j. This result must be unfolded for detector efficiencies and migration 
effects The pair asymmetry A** = 0.071(7) given in Table [3] was computed for dileptonic 

final states with the above cuts, but we expect that it does not change much when computed for 
semileptonic events with the cuts (]IV.2I) . With this proviso, it can be compared with the DO result 
and the numbers agree within experimental and theoretical uncertainties. 

The CDF collaboration measured A" = 0.24 ± 0.14 pS], and the recent result of this experiment 
on the forward-backward asymmetry is A\p^ = 0.193 ± 0.065 ± 0.024 [49]. The CDF results are 
actually unfolded by using SM hypotheses, i.e. they can be compared with the predictions of 
Table [TJ The central values of the CDF results are higher than but still consistent with the above 
SM expectations. The leptonic asymmetries A^ and A^^, which were predicted above at NLOW, 
have not yet been measured to our knowledge. 

Although there is no statistically significant discrepancy between the SM expectations and the 
above cited experimental results, the numbers leave room for speculations about new physics 
contributions. One class of possible new interactions involve axial vector couplings to quarks. At 
the level of uncorrelated tt final states, investigations were made for instance by [15 |I102] . 



IV. 2. Top-spin induced distributions and correlations 

There are a number of angular distributions and correlations with which t, t spin effects, in 
particular tt spin correlations can be analyzed in the ii' and i + j final states. It is well-known 
that the charged lepton from top-decay is, due to the V — A law (which in the meantime is known 
to hold for t ^ 6 to an accuracjill of ~ 10% [39]), the most efficient top-spin analyzer. Thus 
ti spin-correlation effects are most pronounced in i'^i ~ angular correlations in the dilepton final 
states. A correlation of the t or i spin with some axis is best analyzed by a corresponding lepton 
angular distribution in the semileptonic ti modes. 

For qualitative reasoning below we recall here the form of the one-particle inclusive decay density 
matrices integrated over the energies: 

pt-a. =r^^^ + ^^^.^^)^ (IV.8) 
-t^a, ^ rp(^_^^^.^^)^ (IV.9) 

where qi,2 are the directions of flight of ai and a2 in the rest frame of the top and antitop quarks, 
respectively, and the r* are the Pauli matrices. As we incorporate top decay to lowest order in 

^In many SM extensions new-physics contributions to the tbW vertex are quite smah, see e.g. [103] . 
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Figure 12: Left: The rapidity distributions (lIV.Sp . normalized to cr«, for (dashed) and (i~ 
(dotted) for the dileptonic final states at NLOW for the Tevatron for ^ = rrit. Right: The 
differential leptonic charge asymmetry defined in (IIV.4P at NLOW for the Tevatron for n = rrit. 



the electroweak couplings, T^^^ = T^^^ and K2 = ni holds to all orders in as, if 0,2 is the charge- 
conjugate of ai. The spin-analyzing power of leptons is essentially one, K,e+ = Kg- = 0.999, to 
order ag [33]. The spin-analyzing powers of other final state particles/jets are given to order 
in [31]. 

We analyze the following correlations: 
1. The double distributions: 

^ ^'l + Bicos9i + B2Cos92-Ccos9icos92) , (IV.IO) 



a d cos 9 id cos 92 4 

where a denotes the cross section of the respective reaction. The functional form of the 
right-hand side holds only if cuts are not applied. Here 9i (^^2) is the angle between the 
direction of flight of the lepton or jet ji (£'~ or J2) from t {t) decay in the t it) rest frame 
and a reference direction a (b). Below we will choose the axes defined in the last section. 
The coefficient C is a measure of ti spin correlations - it is induced by the correlation ( 1III.5I) . 
In the absence of cuts we have the formula [T3]: 

C = KiK2{ASfk){Sl-h)), (IV.ll) 

which holds for factorizable corrections to all order^. The expectation values on the right- 
hand side of this equation are given in Table HI The coefficients B12 will be discussed in 
items 4 and 5. 



^In fact, the nonfactorizable ©(a^) QCD corrections [35] do not change this result; cf. the comment in jl4] . 
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2. The opening angle distribution [T3] 



--^ = Ul-Dcosv) , (IV.12) 
a a cos Lp 1 

where denotes the angle between the direction of flight of the lepton (or jet ji) and 
of i'~ (or 22)1 defined in the t and t rest frame, respectively. The functional form of the 
right-hand side applies only in the absence of cuts. The coefficient D is induced by a term 
proportional to ■ Sf in the tt production density matrix. If no cuts are applied the formula 

Z} = «:i«:2(^Si-Sf) (IV.13) 

holds. The expectation values on the right-hand side of this equation are given in Table [2l 
Below we comment also on the leptonic opening angle distribution in the laboratory frame. 

3. The dileptonic azimuthal angle correlation, a^^da / dAcf), |511[53]. Here A0 is the difference 
of the azimuthal angles of the and l'~ directions of flight, which can be measured in the 
laboratory frame. This distribution was recently analyzed by [HI] for the case of dileptonic 
tt events at the LHC It was shown by these authors that the A(/) distribution is sensitive to 
SM spin correlations if only events with Mti < 400 GeV are taken into account. The price 
one pays for excluding with this cut a large fraction of the dileptonic events is balanced by 
the fact that the A0 distribution is easier to measure in the dilepton channels than ( IIV.IOI) 
or fHVl2l) . 



4. Nonzero coefficients Bi and B2 (llV.lOP would be due to a correlation St ■ a and Sf ■ b of 
the t and i spin with the respective axis. If polar vectors a, b are chosen then Bi ^ 
would signal parity violation if no cuts were made. For the axes defined in Section IIIII the 
weak-interaction corrections induce non-zero, albeit very small Bi [SniEI]- If cuts are not 
applied one obtains from (IIV.IOI) : 



'^'^ ^:i + 5i,2cos0i,2). (IV.14) 



crd COS 6^1 2 2 

In view of the very small parity-violating effects just mentioned, the SM predicts this dis- 
tribution to be essentially flat |2T] if cuts could be applied that are parity-symmetric with 
respect to the t, respectively i rest frame. However, realistic cuts severely distort (1IV.14I) . 
We compute for the semileptonic ti modes at the LHC: 

^ / X 1 da ^ / X 1 da 

G,+ {z+) = ——, Ge-{z^) = —-^, IV.15 
ai+ dz+ ae- dz_ 

where z± = cos6i± and cos^^+ = a • cos 6*^- = h ■ £~. Below we use for deflniteness the 
helicity basis. One can nevertheless use these distributions for parity and also for CP tests, 
cf. the comments at the end of this section. 
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5. The functional forms of the r.h.s. of Eqs. (1IV.10|) . (]IV.12|) . and (IIV.14P do not hold in 



the presence of kinematic cuts. Kinematic considerations make it clear that the forms 
(]IV.10|) and ( ]IV.14|) will be stronger distorted than (1IV.12|) . Nevertheless, the right-hand 



side of (llV.lOp - or of (1IV.12P - may be used in fits to data for determining the correlation 
coefficients C or D. One may also divide the events in an experimental analysis into spin- 
enhanced (low M^f) and spin-depleted (high M^j) samples and form ratios in order to reduce 
systematic uncertainities |104] . 
Here we use the estimators [H]: 

C = —9 (cos 9i cos 62) , D = —3 (cos ip) , 

5i = 3(cos0i), fi2 = 3(cos02)- (IV.16) 

In the absence of cuts, C = C, D = D, Bi = Bi. Below we compute the expectation values 
(1IV.16P with the above cuts. We drop the hats from now on for ease of notation. In fact, 
we compute the variation of the correlations with Mu- We define it in the case of D by 



dD 

dM^t 



j da cosip5(^sJ{Y,Vfy-Mt?j Qcnt, (IV.17) 



Here 9cut denotes the product of theta functions implementing the cuts and, in case of 
an exclusive computation, of jet functions. For the other angular variables in ( IIV.16|) the 



variations are defined in analogous fashion. We normalize to the respective integrated cross 
section. The computations of (1IV.17I1 and of the other observables concerning their variation 
with respect to Mu made in this paper are theory predictions; i.e. the complications in the 
actual measurement of the invariant mass of the ti pair are not taken into account; cf. the 
comments below. 



Fig. [13] shows dCTacs^^l dMtt and the opening angle distribution flIV.12p for the Tevatron for 



three scales. The resulting integrated correlation coefficients are given in Table HI Comparing 
Cbcam = 0.614(10) with the "no cut" value 0.791(14) shows the significant effects of the cuts on 
this correlation. The distribution dCos/dMu, which is not shown here, has essentially the same 
form as dChcam/ dMu- The correlation coefficient Cofr in the off-diagonal basis is also reduced by 
the cuts to Cos = 0.621(10), cf. Table H Without cuts Cos = 0.798(15). 

The correlations can be understood qualitatively in terms of simple spin arguments. We recall 
this here only for the opening angle distribution, using (1IV.13|) and Born level considerations. The 
ti pair produced by the dominant LO reaction qq — )■ ti is in a ^Si state near threshold. Thus, 
{St ■ St)qq = 1/4 and Dgg = 1/3. The "true" ti intermediate state is of course a mixed state. 
The component from gg fusion, the radiative corrections, and the cuts dilute this correlation to 
D = 0.145(6), cf. TableHl Comparing with the "no cut" value 0.218(2) shows again the significant 
effect of the cuts. 

On may attempt to optimize these correlation coefficients by applying an additional cut on Mu, 
but in view of the not very large number of dileptonic ti events at the Tevatron this seems not 
very useful. In Table H] we give the values of several correlation coefficients for dileptonic events at 
the Tevatron with Mmax = 550 GeV. This value has been chosen here rather arbitrarily; in fact. 
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the correlations do not increase when this cut is tightened to 400 GeV. With even tighter cuts 
one would discard too much dileptonic events. 

Figs. [T^ I15[ and 1161 show dC^ei/dMu, the opening ang le distribution (IIV.12p . and dD/dMu for 
dileptonic final states at the LHC with y/s = 10 and 14 TeV. The distribution dChd/dMu changes 
sign around M^f ~ 700 GeV. This is because for values of M^f not too far above threshold, the 
QCD dynamics tells us that the number of tt pairs produced with like helicities is larger than those 
with unlike helicities, while at large Mu helicity conservation of the strong and weak interactions 
implies that it is the other way around. As the event numbers become small very rapidly and 
because we have normalized to cr«, the distribution dChei/dMu displayed in Figs. [T4l and [TSl is 
essentially zero above ~ 700 GeV. 

As to the qualitative understanding of the opening angle distribution. Figs. [14] and [151 Because 
the production of ti pairs at the LHC is dominated by gluon fusion, the pair is, to LO QCD, 
in a ^Sq state at threshold. Hence (St ■ Sf)^^ = —3/4 and Dgg = —1. Away from threshold, 
the pair is actually in a mixed state already at Born level. This fact, the component from qq 
annihilation, and the radiative corrections deplete this correlation to D = —0.252(5) (10 TeV), 
respectively D = —0.240(10) (14 TeV), cf. Table [H The effect of cuts on D is less drastic than 
in the case of the Tevatron, as seen by comparing with the "no cut" values —0.233(4) (10 TeV) 
and —0.236(3) (14 TeV). On the other hand, the helicity correlation is depleted by the cuts, as 
shown by comparing the respective numbers in Table [H with Chei(no cut) = 0.326(2) (10 TeV) and 
0.328(3) (14 TeV). 

The helicity and opening angle correlations at the LHC can be enhanced significantly by taking 
into account only events with a tt invariant mass below a certain value Mmax- Table [H shows that 
for Minax = 550 GeV the correlation coefficients D and Chci can be enhanced by about 43% and 
50%, respectively. We have not optimized this cut. 

The beam, off-diagonal, and helicity correlation and the opening angle distribution were computed 
before to NLO QCD for dileptonic, semileptonic, and non-leptonic final states in [14J, without cuts, 
and with a different set of PDF (CTEQ6.1M and MRST2003) and a slightly different value of 
the top-quark mass. The numbers given in that paper can be compared with the results given 
here, using Table [2l the formulae flIV.ll|) and (IIV.13|) . and the spin-analyzing powers k given 
in [H]. There are two effects: the PDF set CTEQ6.6M and the smaller value of the top-quark 
mass that we use in the present paper enhances the correlation coefficients slightly as compared 
to those of [H]. On the other hand the weak-interaction corrections lead to a slight reduction of 
the correlations at the LHC (cf. Sect. while their infiuence on the Tevatron correlations is 
smaller. In effect this leads to slightly enhanced numbers for the Tevatron, about 2%, while the 
numbers for the LHC given in Table [2] are essentially the same as those in [H] . 
The authors of jHH] computed the opening angle distribution as defined above in the MC@NLO 
framework - whith parton showering, but no cuts were applied - for dileptonic final states at the 
Tevatron and the LHC (14 TeV) and found agreement with the results given in [T5|. This shows 
that this leptonic correlation is essentially unaffected by the parton shower. This distribution was 
recently also evaluated in [91] for dileptonic final states at the Tevatron and the LHC (10 TeV), 
to NLO QCD with cuts. We have computed these distributions also with the selection criteria 
of ^9]J and agree with their results. 
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Figure 13: Distributions at NLOW for ^^ final states at the Tevatron: dCi,eam/dMtt (left panel) 
and the opening angle distribution (1IV.12P (right panel). The solid, dashed, and dotted lines 
correspond to ^ = rrit, mt/2, and 2mt, respectively. 

The dileptonic angular correlations considered above require the reconstruction of the t and t rest 
frames, which is a highly non-trivial experimental task. There are enough kinematic constraints 
but, as is well-known, for dileptonic events a four-fold ambiguity arises, which may increase to an 
eightfold one if the b and b jets are incorrectly assigned. Monte-Carlo based methods have been 
devised (see, e.g., |105j ) in order to reduce or partly resolve these ambiguities. 
The beam-basis and off-diagonal correlations were recently measured for dileptonic events at the 
Tevatron by the DO and CDF collaborations, respectively, with the results 

Cbeam = -0.17l°|^ [65], /toff = 0.32+°;^ m- (IV.18) 

In our sign convention Kos = Cos- These results agree, within the still large experimental uncer- 
tainties, with our predictions for the dilepton channels at the Tevatron given Table HI The CDF 
collaboration measured also the helicity correlation for lepton + jets events. They obtained for 
the respective spin correlation parameter at the ti level the value n^ei = 0.60 ± 0.50 ± 0.16 [67] . 
In our sign convention for the double spin asymmetries [H] this is equivalent to — Chei and should 
be compared with the value — (9hei = +0.368(10) for the Tevatron given in Table [2l The above 
acceptance cuts deplete this value to — Chci = +0.30(6), cf. Table IH 

A simulation study ^106J within the ATLAS collaboration concluded that D and Chei can be mea- 
sured at the LHC with a relative uncertainty of 6D ^ 5% and SC^ei — 7%, respectively, using 
dileptonic and semileptonic ti events. 
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Mjj (GeV) cos(p„ 

Figure 14: Distributions at NLOW for ii final states at tfie LHC {y/s = 10 TeV): dC^^i/dMu (left 
panel) and the opening angle distribution (IIV.12P (right panel). The solid, dashed, and dotted 
lines correspond to = rrit, mt/2, and 2mt, respectively. 




M„ (GeV) cos(p„ 

Figure 15: Distributions at NLOW for ££ final states at the LHC {y/s = 14 TeV): dChei/dMa (left 
panel) and the opening angle distribution (IIV.12P (right panel). The solid, dashed, and dotted 
hues correspond to /i = m^, mt/2, and 2mt, respectively. 
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Figure 16: The distribution dD/dMa at NLOW for U final states at the LHC at = 10 TeV 
(left panel) and at 14 TeV (right panel). The solid, dashed, and dotted lines correspond to 
= mt,mt/2, and 2mt, respectively. 



Table 4: Results for several observables at NLOW for dileptonic final states at the Tevatron and 
LHC, with cuts as specified in ( IIV.II) . The correlation coefficients are computed also with a cut 
Mtt < Mmax- Here we have chosen M^ax = 550 GeV. 





Tevatron 


LHC (10 TeV) 


LHC (14 TeV) 




mt/2 


mt 


2mt 


mt/2 


mt 


2mt 


nit/2 


m 


2m 


CTu (pb) 


0.043 


0.042 


0.038 


2.31 


2.03 


1.76 


5.00 


4.38 


3.82 


D 


0.139 


0.145 


0.151 


-0.257 


-0.252 


-0.257 


-0.240 


-0.247 


-0.230 




0.125 


0.132 


0.138 


-0.344 


-0.340 


-0.347 


-0.340 


-0.353 


-0.338 


Chel 


-0.294 


-0.299 


-0.306 


0.249 


0.247 


0.252 


0.225 


0.237 


0.229 


CheKAfmax) 
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Next we analyze two leptonic angular correlations defined in the laboratory frame. The distribu- 
tion of the opening angle between the lepton momenta in the laboratory frame, da /d cos ip^ , 
which was recently considered at NLO QCD in [91j, is almost insensitive to tt spin correlations [96 . 
This is demonstrated in Fig. [T7] (left panel) where this distribution is shown for the LHC with|j 
and without tt spin correlations. The correlation of the boosted leptons in the laboratory frame 
traces the correlation of the t and t spins inefficiently. On the other hand, the distribution (IIV.12P 
discriminates by constructioij^ between spin correlations "switched-on" and "switched-off " , cf. 
Fig. [TT] (right panel). 

As shown recently [64] by a tree-level analysis, the dileptonic azimuthal angle correlation at the 
LHC, a~^da/dA(j), does discriminate between correlated and uncorrelated ti events if dileptonic 
final states with Mu < 400 GeV are taken into account only. Here A0 is the difference of the 
azimuthal angles of and i'~ in the laboratory frame. This distribution is plotted in Fig. [T8]at LO 
and NLOW for correlated and uncorrelated dileptonic tt events. In [61] only the LO distribution 
is given with which we agree. For the cross section ratio r = au{Mu < 400 GeV)/aee we get with 
the above cuts r ~ 18.6% at NLO QCD with /i = nit. Assuming an integrated luminosity of 10 
fb~^ at the LHC (14 TeV), this corresponds to ~ 32000 dilepton events {i = e,ij). 
The shapes of the correlated and uncorrelated distribution depend very sensitively on the cut on 
Mff. Thus, a critical (experimental) issue will be how well the theory-cut on can be realized in 
the selection of dileptonic events at the LHC, given the ambiguities that arise from the presence 
of the two neutrinos in the final state. This complication was analyzed in [M]. In addition, the 
azimuthal angle correlation looses its discriminating power the more the cut on is relaxed 
to values above Mmax = 400 GeV. Fig. [19] shows this distribution at NLOW for correlated and 
uncorrelated dileptonic tt events (with cuts as given in (IIV.1|) ) if no cut on Mu is applied. The plot 



compares well with Fig. 1 of [88] computed with the MC@NLO code. (The above comments on the 
differences between [HS] and our approach in the computation of the openening angle distribution 
apply also here.) 

The azimuthal angle correlation is clearly easier and more precisely measurable (assuming that the 
problem of finding the best option for selecting dileptonic events will be solved) than the helicity 
correlation or the opening ang le distribution ( IIV.12p . It can be used to check the SM prediction 
on tt spin correlations or may be employed in the search for new physics effects in the low energy 
tail of the Mu spectrum. If one one wants to keep spin correlation observables in the tool kit, 
especially for new physics searches in the upper part of the Mu spectrum, where new resonances 
or continuum effects may show up, one should use observables like the helicity correlation or the 
opening angle distribution. 

The beam, off-diagonal, and helicity correlation and the opening angle distribution (1IV.12P can 
also be measured for semileptonic tt events at the Tevatron and at the LHC. In this case the t and t 
rest frames are easier and more efficiently reconstructible than for dileptonic final states. Moreover 
the number of i+ jets events is about six times larger than the number of ii' events. In the case of 
a hadronically decaying top quark the least energetic non-6 jet j< among the top-decay products 



^For comparison we have computed this distribution also with the cuts used in [91 and find agreement with 
Fig. 4 of this work. 
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This is of course also the case for the beam, off-diagonal, and helicity correlations. 
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Figure 17: Left: The opening angle distribution a^^da / dcosip^ in the laboratory frame. Right: 
The opening angle (IIV.12P defined with respect to the tt ZMF. Both are computed at NLOW with 
fi = rrit for ii final states at the LHC {^/s = 14 TeV). The solid lines are the SM predictions, the 
dashed lines are the distributions which result when the tt spin correlations are switched off. 

may serve as analyzer of the top spin. It was shown [35] that j< has the best spin-analyzing power 
for this decay mode (barring the extremely difficult task of tagging the flavor of the d-tjpe jet from 
W decay); Kj^ ~ 0.47 to order depending on the jet algorithm. Thus the angular correlations 
for the i+ jets events will be reduced roughly by this factor as compared with the respective cor- 
relations for ii' final states. Detailed predictions with acceptance cuts will be given elsewhere |107] . 

Finally we consider, for the the i'^+ jets modes, the one-particle distributions G£±{z±) defined in 
(IIV.ISP where z± = cos 9i± denote the cosines of the lepton helicity angles, i.e., 9^+ {Oi-) is the 
angle between the i^ {i~) direction of flight in the t it) rest frame and the t it) direction in the 
tt ZMF. 

If no cuts are applied, non-trivial, i.e. non-fiat distributions would arise if the ensemble of (anti)top 
quarks were longitudinally polarized - a small effect in the SM 1201121]. However, the presence of 
cuts severely distort these distributions in the vicinity of cos9i± = — 1. This region corresponds 
to leptons that are emitted into the hemisphere opposite to the (anti)top direction of flight in 
the ti ZMF. These leptons are therefore less energetic on average and thus stronger affected by 
cuts than those in the remaining region. This is shown in Fig. [20] for and the corresponding 
differential distribution dBi/dMti at the LHC. The SM value of Bi is given, for the above cuts, 
in Table m To NLO in the gauge couplings, the distributions Gg~ and dB2/dMu are identical to 
G£+ and dBi/dMu, respectively, which we checked with our code. This is because 

G,+ (^+) = G,-(^_) (IV.19) 
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Figure 18: The azimuthal angle distribution a^^{M„-^ax = 4:00GcY)da/dA(f) in the laboratory 
frame for ii final states at the LHC (14 TeV) with and without ti spin correlations, with the 
additional cut Mu < M^nax- Sohd and dashed hues are LO and NLOW, respectively, for /j, — nit. 
Right: The distribution at NLOW for ji — rrit (solid), mt/2 (dashed), and 2mt (dotted). 




Figure 19: Left: The azimuthal angle distribution a~^da / d^cf) in the laboratory frame for 11 final 
states at the LHC (14 TeV) at NLOW for /j, = rrit and no cut on M^f. The solid line is with and 
the dashed line without tt spin correlations. 
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Figure 20: The distribution a^^da / dcos6i+ (left panel) and the corresponding differential distri- 
bution dBi/dMtt (right panel) for £++ jets at the LHC (14 TeV). Here 6^+ is the polar angle of 
in the t rest frame with respect to the t direction in the ti ZMF. The NLOW results are for 
fi = nit (sohd), mt/2 (dashed), and 2mt (dotted). 

holds in the SM to a high degree of accuracy if the selection cuts/criteria are CP images of each 
other. Although the pp state is not a CP eigenstate, one can show that (]1V.19|) holds as long as 



gg and qq initiated ti production is governed by CP-invariant interactions [112. . 
One may use these distributions for checks of whether or not the collected ensemble of (anti)top 
quarks has a longitudinal polarization, which would be due to some new parity-violating inter- 
action that affects ti production^. That is, one may investigate whether or not the measured 
distributions differ from the SM results. Fig. |20l 

If one defines cos6e+ = a ■ with respect to an axial vector a, e.g. the vector normal to the 
2 — )■ 2 scattering plane (and likewise for one may check for a normal polarization of top 
quarks generated by absorptive parts in the scattering amplitude. At NLO QCD this polarization 
amounts to a few percent [58|l59]. 

We close this section with a few remarks on some distributions/correlations suitable for CP sym- 
metry tests respectively searches for non-standard CP violation (cf., e.g., |109] for a review). This 
can be made with dileptonic and semileptonic ti events at the Tevatron and especially at the LHC, 
once sufficiently large data samples will have been collected. Although the (expectation values of 
the) observables below cannot be classified a priori with respect to CP in the case of pp collisions, 
they may nevertheless be used at the LHC. As shown in |110yillj by a general kinematic analysis, 
P- and CP- violating interactions that affect hadronic ti production induce, in general, at the level 



Related parity-violating polarization observables and their size in some SM extensions were investigated in |108| . 
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of the tt intermediate states two types of spin correlations/asymmetries: 

1. A difference in the longitudinal polarization of t and t quarks. The most useful asymmetry 
corresponds to k ■ (S^ — S^) where k is the direction of the t quark in the tt ZMF. In the 
dileptonic or semileptonic tt events these correlations lead to a difference in the distributions 
Gi±. One may consider the ratio 

As shown above, A^p = to NLOW in the SM. 

A non-zero A^p requires, apart from P and CP violation, also an absorptive part in the 
scattering matrix. A non-zero polarization asymmetry k ■ {St — St) amounts to a difference 
N{tLtL) — Nitniji) 7^ of the number of ti with negative and positive helicity |113j . In the 
case of a CP-violating Higgs sector, in particular if heavy spin-zero states with undefined 
CP parity exist that are resonantly produced in the tt channel, then suitable asymmetries 
evaluated in suitable bins around the position of the resonance [HS] may be as large as a few 
percent [112j . 

In the context of leptonic charge asymmetries discussed in Sect. llV.ll it is natural to consider 
or its integrated version, 

/ dyN,+ {y)- J dyN,-{y) 
"^'^^ = '7dyN,,iy)+JdyN,-iy) ' ^^^"''^ 

y>0 y<0 

The asymmetries (]IV.21I) and (]IV.22I) correspond to a top-spin asymmetry with respect to 
the beam direction. They are rather insensitive to P- and CP-violating Higgs-boson effects. 
Larger effects would result from s-channel exchanges of new heavy J 7^ resonances with 
CP-violating couplings to top quarks. A possible difference in the transverse energies of 
and i~ , which derives also from these single spin asymmetries, is below the percent level for 
non-resonant CP- violating exchanges of Higgs bosons |113j or of squarks/gluinos [114]. 

2. P- and CP-odd/T-odd ti spin correlations. The most useful one is proportional to k-(St x Sf). 
It can be induced already at tree level. In dileptonic events it would show up in correlations 
proportional to Ocp = k ■ x £~), and likewise for semileptonic events. When the unit 
momenta and k are taken in the t (t) rest frame and in the tt ZMF, respectively, 
then Ocp corresponds, in the absence of cuts, to a term proportional to sin A0*, where A0* 
denotes the difference of the i~ azimuthal angles in the plane orthogonal to k. As the 
construction of these frames with reasonably small experimental error will be difficult, one 
may use pseudo rest-frames in analogy to what was shown in |115j . 

In the case of resonant s-channel exchange of a Higgs particle with scalar and pseudoscalar 
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Yukawa couplings, suitable asymmetries Acp oc N{Ocp > 0) — N{Ocp < 0) may be 
as large as a few percent when evaluated in appropriate bins around the position of the 
resonance [9611112] . and do not receive SM contributions at the few per mille level. 

The above dileptonic azimuthal angle correlation in the laboratory frame, a^^da / d/\(j), is 
symmetric with respect to A0 — t- 27r — A0 in the case of CP invariance. CP-violating 
interactions would, in general, modulate it by a term proportional to sin A0 (in the absence 
of cuts); i.e. one should measure it in the whole range < A0 < 27r. In the terminology 
of [llOj this term corresponds to a correlation of S( x Sf with the beam direction in the 
laboratory frame, i.e., to a correlation observable pl ■ (£j x £2). As far as possible Higgs 
sector CP violation is concerned, this correlation is rather insensitive to resonant terms if the 
mass of the resonance is significantly above the tt production threshold. The non-resonant 
Higgs-boson exchange terms lead to rather small CP-violating effects [ITIj . 

V. Conclusions 

We have investigated a number of observables that are and will be instrumental in the exploration 
of tt production and decay at the Tevatron and the LHC. For this analysis we made a computer 
program that incorporates besides the NLO QCD corrections to ti production and decay also 
mixed weak-QCD corrections to the production amplitudes, and that allows for studies of corre- 
lated versus uncorrelated ti events. 

We assessed, for a few distributions and correlations at the LHC, the relative size of the weak- 
interaction contributions as compared to the NLO QCD results. In the px or range where most 
of the tt events at the LHC are located, these contributions deplete the NLO QCD results some- 
what, for instance, the px and M^i distribution by about —2% for pt ~ 400 GeV and M^^ ~ 1.2 
TeV. For larger pt, M^i the depletion grows to several percent. As to the charge asymmetries A 
and A** computed at the level of intermediate tt: our results agree with [IS]; yet we find that the 
weak-interaction contributions are slightly smaller than those determined in [13] , due to additional 
terms that were previously not taken into account. 

Our main physics results are the investigations at NLOW in the gauge couplings of the asymme- 
tries and angular correlations for dileptonic events which were presented in Sect. IIVI First, we 
considered charge asymmetries at the Tevatron. The pair- asymmetry A** decreases by about 10% 
if acceptance cuts are applied. Our result, A** = 0.071(7) (where the given uncertainty due to 
scale variations certainly underestimates the true theoretical error), may be compared with the 
respective DO measurement [17]. Both results agree within the rather large uncertainties. We 
computed also a leptonic charge asymmetry and a pair-asymmetry A^^ for the Tevatron, which 
are related to A and A**, respectively. As in the case of A**, we determined these asymmetries 
for correlated and uncorrelated ti events, in order to study the effect of spin correlations on these 
distributions. We showed that ti spin correlations affect the pair asymmetry A^^ by about 7% for 
the acceptance cuts used in Sect. [IVl These leptonic asymmetries have, to our knowledge, not 
yet been measured, probably because of limited numbers of (di) leptonic events. As these leptonic 
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asymmetries should eventually be measurable more precisely than A or A**, their experimental 
determination might provide a more conclusive comparison with SM results. Our NLOW predic- 
tions are given in Table [3l One may speculate that top-spin effects are more pronounced in more 
exclusive charge or forward-backward asymmetries [S], but this remains to be investigated. 

Furthermore, we determined several dileptonic angular correlations, which reflect tt spin correla- 
tions with respect to different spin bases, namely the beam, off-diagonal and helicity correlation, 
and an opening angle distribution defined in a specific way [H], when selection cuts are applied. 
Our NLOW predictions for the beam, off-diagonal, and helicity correlation for the Tevatron agree 
with recent measurements by the CDF and DO experiments; albeit the experimental uncertainties 
are still large, mainly because of limited statistics. 

We have also made predictions for estimators (which have 1-dimensional distributions) of these 
correlations as functions of Mu. These estimators may prove useful for the Tevatron and also in 
the early rounds of LHC data-analyses, where the event numbers will likely not be abundant. 
We have also considered two angular correlations defined in the laboratory frame for dileptonic 
final states at the LHC, and we have computed them to NLOW for correlated and uncorrelated 
ti events. Our results confirm the findings of [61], namely that the azimuthal angle correlation is 
a sensitive variable for exploring tt spin correlations at the LHC in the low-energy tail of the 
spectrum (Mjf < 400 GeV). It should be recalled that the helicity correlation and the opening 
angle distribution as defined in [13] are sensitive to spin correlations also when applied to event 
samples with M^f^^ > 400 GeV. An crucial experimental issue in the precise measurement of both 
the azimuthal angle correlation and the latter observables is the kinematic reconstruction of the 
dileptonic events. 

Moreover, we have briefly discussed a few P- and CP-odd observables that should be useful in 
the search for P- and/or CP- violating interactions in dileptonic and semileptonic tt events. Our 
NLOW results on distributions and correlations indicate that the above observables do not receive 
SM contributions at the level of a few per mille or below. This issue deserves further, more detailed 
studies. 

It is also worth recalling here that, in case the measurements of the tt spin correlations at the LHC 
should eventually match the SM predictions - for instance predictions that involve the low-energy 
tail of the spectrum - and would reach a reasonable level of precision, one may use these 
correlations also for the further exploration of the parton content of the proton. They are, in fact, 
quite sensitive to the proton's relative quark and gluon content, because the contributions from 
gg and qq initiated tt production have opposite signs [T^IT^. 

The above investigations and comparisons with existing experimental results show that the pre- 
cision, with which some distributions, especially angular correlations can be predicted, will likely 
not be matched quickly by experiments in the near future. Nevertheless, in view of the poten- 
tial importance of observables such as the leptonic charge asymmetries and spin correlations in 
exploring the interactions of top quarks, further detailed investigations on expected measurement 
uncertainties at the LHC would certainly not be futile exercises. 
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Future extensions of our work may include the implementation of threshold resummations into 
our code and, as to more phenomenological issues, an analysis of charge asymmetries and angular 
correlations at NLOW also for lepton + jets final states, and a more systematic error analysis by 
taking into account also other available PDF sets |116[lll7j including the PDF uncertainties. 
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